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Abstract 


This  manuscript  is  the  first  in  a  sequence  of  documents  that  describe  the  development  and 
verification  of  the  various  components  of  a  sea  surface  infrared  radiance  simulator.  We  have 
developed  an  approach  to  modeling  the  sea  surface  radiance  in  the  infrared  to  include  both 
a  coherent  ship  wake  structure  and  an  averaged  background  component.  In  this  manuscript 
we  discuss  the  component  of  the  simulator  for  the  ambient  sea  surface.  The  contribution 
from  the  background  sea  surface  is  determined  by  its  roughness  and  temperature  given  a 
defined  incident  sky  radiation.  The  roughness  is  largely  determined  by  the  wind  speed  at 
the  sea  surface  and  is  described  as  a  probability  distribution  of  slopes  of  planar  facets  in  the 
model.  Then  given  the  reflectivity  and  emissivity  of  the  sea  surface  and  their  inherent  de¬ 
pendencies  on  angle  to  the  observer,  temperature,  and  roughness  we  calculate  the  infrared 
radiance  from  the  sea  surface.  We  have  illustrated  the  range  of  behaviour  of  the  sea  surface 
radiance  for  simple  sky  radiances  at  grazing  and  non-grazing  angles  as  a  function  of  sea 
roughness  and  temperature. 


Resume 


Ce  rapport  est  le  premier  d’une  serie  de  documents  qui  decrivent  le  developpement  et  la 
verification  de  plusieurs  composantes  d’un  simulateur  de  la  radiance  de  la  surface  de  l’eau 
dans  l’infrarouge.  Nous  avons  developpe  une  approche  pour  modeliser  la  radiance  de  la  sur¬ 
face  de  l’eau  dans  l’infrarouge  incluant,  la  fois,  la  contribution  moyenne  de  Tarriere-plan 
maritime  ainsi  qu’une  structure  coherente  de  sillages  de  navires.  Dans  le  present  document, 
nous  traitons  la  composante  du  simulateur  qui  conceme  la  contribution  de  la  radiance  de 
la  surface  de  l’eau  de  Tarriere-plan  maritime.  Cette  radiance  emise  par  la  surface  de  l’eau 
est  determinee  principalement  par  sa  temperature  et  sa  rugosite.  La  rugosite  de  la  sur¬ 
face  de  l’eau  depend  largement  de  la  vitesse  du  vent  a  sa  surface.  Elle  est  decrite  par  une 
fonction  de  densite  de  distribution  de  probability  des  pentes  des  differentes  facettes  pla- 
naires  qui  represented  la  surface  de  l’eau  dans  le  model.  Ainsi,  connaissant  l’emissivite 
et  la  reflectivite  de  la  surface  de  l’eau,  leur  dependance  avec  Tangle  de  l’observateur,  la 
temperature  de  l’eau  et  sa  rugosite,  nous  avons  calcule  la  radiance  emise  dans  l’infrarouge 
de  la  surface  de  l’eau  et  celle  qui  y  est  reflechie  d’un  ciel  donnee.  Nous  avons  illustre  la 
dependance  du  comportement  de  la  radiance  de  la  surface  de  l’eau  pour  differents  angles 
d’observateur  en  fonction  de  la  rugosite  et  de  la  temperature  de  la  surface  de  l’eau  en  utili- 
sant  un  modele  simplifie  de  la  radiance  du  ciel. 
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Executive  summary 


Sea  Surface  Infrared  Radiance  Simulator:  Part  1 : 
Roughness  and  Temperature  models  of  the  Sea 
Surface  Radiance 

Vivian  Issa,  Zahir  A.  Daya;  DRDC  Atlantic  TM  2010-280;  Defence  R&D  Canada  - 

Atlantic;  December  2010. 

Background:  Images,  and  the  data  extractable  from  them,  are  an  indispensible  component 
of  a  modern  defence  analysis  capability.  In  the  maritime  domain,  the  sea  surface  which  is 
seldom  stationary,  has  a  broad  range  of  appearances  in  the  visible  and  infrared  optical 
bands.  The  optical  appearance  of  the  sea  surface  is  influenced  by  wind,  sky  and  clouds 
so  much  so  that  the  basic  measureables  such  as  brightness  or  contrast  vary  significantly 
with  prevailing  and  changing  conditions.  This  sensitive  dependence  may  be  exacerbated 
in  the  infrared  where  the  sea  surface  and  the  atmosphere  are  themselves  thermal  sources. 
There  are  several  approaches  to  modeling  the  sea  surface  and  its  optical  appearance  visually 
and  in  the  infrared,  however,  given  the  perpetual  motion  and  inherent  stochasticity  of  the 
sea  surface,  the  measureables  are  invariably  temporal  or  spatial  averages.  Yet  there  are 
often  deterministic  or  long-lived  coherent  features  on  the  sea  surface  such  as  ship  wakes, 
swells  and  gravity  waves.  Ship  wakes,  in  particular,  are  important  features  in  airborne  and 
satellite-based  surveillance  and  are  speculated  to  be  a  defence  vulnerability. 

Principal  results:  We  have  developed  an  approach  to  modeling  the  sea  surface  radiance 
in  the  infrared  to  include  both  a  coherent  ship  wake  structure  and  an  averaged  background 
component.  The  ship  wake  structure  will  consist,  geometrically,  of  a  Kelvin  gravity  wave 
field,  a  white  water  wake  and  a  turbulent  trailing  wake.  The  background  component  derives 
from  the  probabilistic  description  of  slopes  of  the  sea  surface  due  to  the  wind.  Combining 
the  sky  and  atmospheric  radiation  that  is  reflected  by  these  features  on  the  sea  surface  with 
the  emitted  radiation  from  the  water  and  the  radiation  along  the  path,  we  can  model  the 
total  radiance  received  by  the  observer.  In  brief,  our  sea  surface  simulator  will  consist  of 
modular  components  for  the  ambient  sea  surface,  for  the  ship  wake,  for  the  sky  radiance 
and  for  the  extended  source  to  finite  receiver  geometry.  At  the  present  time,  we  discuss  the 
component  of  the  simulator  for  the  ambient  sea  surface. 

Significance:  We  have  built  a  tool  to  simulate  the  reflected  and  emitted  radiance  from  the 
sea  surface  knowing  its  temperature  and  geometry.  This  tool  is  a  fundamental  step  for  the 
study  of  the  ships  wakes  radiance.  Wakes  are  important  because  they  are  signatures  rich 
in  information  and  they  persist  on  the  surface  of  the  water  over  several  kilometers,  and  for 
several  minutes.  This  importance  is  even  more  pronounced  in  the  infrared  because  the  sea 
surface  and  sky  are  thermal  sources. 
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Future  work:  This  manuscript  is  the  first  in  a  sequence  of  documents  that  describe  the 
development  and  verification  of  the  various  components  of  our  sea  surface  infrared  radiance 
simulator.  Subsequent  manuscripts  will  describe  the  discretization  and  mapping  of  the  sea 
surface  to  the  observer’s  receiver  optics,  the  geometrical  and  infrared  models  for  the  ship 
wake  and  the  performance  of  the  simulator. 
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Sea  Surface  Infrared  Radiance  Simulator:  Part  1 : 
Roughness  and  Temperature  models  of  the  Sea 
Surface  Radiance 

Vivian  Issa,  Zahir  A.  Daya ;  DRDC  Atlantic  TM  2010-280 ;  R  &  D  pour  la  defense 

Canada  -  Atlantique ;  decembre  2010. 

Contexte  :  Les  images  ainsi  que  les  donnees  que  l’on  peut  extraire  d’elles  sont  des  elements 
indispensables  pour  les  analyses  de  la  defense  moderne.  Dans  le  domaine  maritime,  la  sur¬ 
face  de  l’eau,  qui  est  rarement  stationnaire,  a  une  grande  gamine  d’apparences  dans  les 
bandes  optiques  du  visible  et  de  l’infrarouge.  Le  vent,  le  ciel  et  les  nuages  influencent  gra- 
dement  l’apparence  optique  de  la  surface  de  l’eau  d’ou  la  grande  dependance  des  quantites 
mesurables,  telles  que  la  brillance  et  le  contraste,  avec  ces  parametres.  Cette  dependance 
est  accentuee  dans  l’infrarouge  ou  le  ciel  et  la  surface  de  l’eau  sont  des  sources  thermiques. 
II  existe  plusieurs  approches  pour  modeliser  la  surface  de  l’eau  et  son  apparence  optique 
dans  le  visible  et  dans  l’infrarouge.  Etant  donnee  le  mouvement  perpetuel  et  stochastique 
de  la  surface  de  l’eau,  ces  variables  sont  moyennement  constant  temporellement  et  spa- 
tialement.  Toutefois,  la  surface  de  l’eau  peut  avoir  des  traits  coherentes  deterministes  ou 
d’une  duree  plus  longue  tels  que  les  sillages  des  navires,  la  houle  et  les  vagues  de  gravites. 
Particulierement,  les  sillages  des  bateaux,  sont  des  traits  importants  pour  la  surveillance 
satellite  et  aerienne,  ce  qui  est  d’une  grande  importance  pour  les  applications  de  la  defense. 

Resultats  principaux  :  Nous  avons  developpe  une  approche  pour  modeliser  la  radiance 
de  la  surface  de  l’eau  dans  1’ infrarouge  pour  inclure  a  la  fois  l’arriere  plan  maritime  et 
les  sillages  des  navires.  La  structure  du  sillage  est  geometriquement  constitute  d’une  onde 
de  gravite  de  Kelvin,  d’un  sillage  d’ecume  (zone  d’eau  blanche)  et  du  sillage  turbulent. 
L’arriere  plan  maritime  est  affecte  par  la  vitesse  du  vent.  Cette  derniere  est  un  parametre 
important  dans  la  fonction  de  distribution  de  densite  de  probability  des  pentes  des  facettes 
de  la  surface  de  l’eau.  La  radiance  du  ciel  et  de  l’atmosphere  reflechie  sur  ces  facettes  ainsi 
que  la  radiance  intrinseque  emise  par  ces  memes  facettes  sont  les  deux  composantes  de  la 
radiance  regue  par  un  observateur.  Bref,  pour  une  geometrie  specifique  de  la  source  et  de 
l’observateur,  notre  simulateur  de  la  radiance  de  la  surface  de  l’eau  va  contenir  plusieurs 
composantes,  soient  la  radiance  de  la  surface  de  l’eau  ambiante,  celle  de  sillages  de  navires 
ainsi  que  celle  de  la  radiance  du  ciel.  Presentement,  nous  traitons  la  composante  de  la 
radiance  de  la  surface  de  l’eau  ambiante. 

Importance  :  Nous  avons  un  outil  pour  simuler  la  radiance  reflechie  et  emise  par  la  surface 
de  l’eau  sachant  sa  temperature  et  sa  geometrie.  Cet  outil  est  fondamentale  pour  l’etude  de 
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sillages  des  navires.  Les  sillages  des  navires  sont  des  signatures  importantes  puisqu’elles 
sont  riches  en  informations  et  persistent  sur  la  surface  de  l’eau  sur  plusieurs  kilometres 
et  pour  plusieurs  minutes.  Cette  importance  est  encore  plus  prononcee  dans  l’infrarouge 
puisque  le  ciel  et  la  surface  de  l’eau  sont  des  sources  thermiques. 

Recherches  futures  :  Ce  document  est  le  premier  d’une  serie  de  documents  qui  decrivent 
le  developpement  et  la  verification  de  plusieurs  composantes  d’un  simulateur  infrarouge  de 
la  radiance  de  la  surface  de  l’eau.  Des  documents  ulterieurs  suivront  pour  decrire  le  modele 
de  discretisation  de  l’espace  geometrique,  le  modele  infrarouge  de  sillages  des  navires  ainsi 
que  la  performance  du  simulateur. 
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1  Introduction 


Images,  and  the  data  extractable  from  them,  are  an  indispensible  component  of  a  modern 
defence  analysis  capability.  In  the  maritime  domain,  the  sea  surface  which  is  seldom  sta¬ 
tionary,  has  a  broad  range  of  appearances  in  the  visible  and  infrared  optical  bands.  The 
optical  appearance  of  the  sea  surface  is  influenced  by  wind,  sky  and  clouds  so  much  so  that 
the  basic  measureables  such  as  brightness  or  contrast  vary  significantly  with  prevailing  and 
changing  conditions.  This  sensitive  dependence  may  be  exacerbated  in  the  infrared  where 
the  sea  surface  and  the  atmosphere  are  themselves  thermal  sources. 

The  sea  surface  is  complex:  the  surface  is  not  constant,  it  is  not  always  mathematically 
smooth  having  cusps  and  sigularities  when  it  is  rough,  it  extends  to  the  limits  of  visibility, 
its  appearance  depends  on  the  position  of  the  observer  (the  viewing  angle  from  the  average 
normal  to  the  sea  surface)  reflecting  different  parts  of  the  sky  which  can  be  clear  or  cloudy 
etc.  There  are  several  approaches  to  modeling  the  sea  surface  and  its  optical  appearance 
visually  and  in  the  infrared  [1,  2,  3,  4,  5],  however,  given  the  perpetual  motion  and  inherent 
stochasticity  of  the  sea  surface,  the  measureables  are  invariably  temporal  or  spatial  aver¬ 
ages.  Yet  there  are  often  deterministic  or  long-lived  coherent  features  on  the  sea  surface 
such  as  ship  wakes,  swells  and  gravity  waves.  Ship  wakes,  in  particular,  are  important 
features  in  airborne  and  satellite-based  surveillance  (see  e.g.  [6])  and  are  speculated  to  be 
a  defence  vulnerability. 

Our  approach  to  modeling  the  sea  surface  radiance  in  the  infrared  is  designed  to  include 
both  a  coherent  ship  wake  structure  and  an  averaged  background  component.  The  ship 
wake  structure  will  consist,  geometrically,  of  a  Kelvin  gravity  wave  field,  a  white  water 
wake  and  a  turbulent  trailing  wake.  The  background  component  derives  from  the  proba¬ 
bilistic  description  of  slopes  of  the  sea  surface  due  to  the  wind.  Combining  the  sky  and 
atmospheric  radiation  that  is  reflected  by  these  features  on  the  sea  surface  with  the  emitted 
radiation  from  the  water  and  the  radiation  along  the  path,  we  can  model  the  total  radiance 
received  by  the  observer.  Whereas,  this  general  ‘summation  of  different  parts’  approach 
appears  simple  in  principle,  its  implementation  is  rife  with  intricacies,  some  of  which  are 
discussed  in  Section  3.  In  brief,  our  sea  surface  simulator  will  consist  of  modular  com¬ 
ponents  for  the  ambient  sea  surface,  for  the  ship  wake,  for  the  sky  radiance  and  for  the 
extended  source  to  finite  receiver  geometry. 

This  manuscript  is  the  first  in  a  sequence  of  documents  that  describe  the  development  and 
verification  of  the  various  components  of  our  sea  surface  infrared  radiance  simulator.  In  this 
manuscript,  we  briefly  describe  the  overall  framework  (Section  3)  and  the  module  for  the 
infrared  radiance  from  the  ambient  or  background  sea  surface.  This  module  accounts  for 
the  reflected  infrared  radiance  from  the  sky  (Section  4)  and  for  thermally  emitted  infrared 
radiance  from  the  water  (Section  5).  Results  of  the  numerically  computed  infrared  radiance 
as  a  function  of  sea  surface  roughness,  temperature  and  viewing  angle  are  discussed  in 
Section  6  which  is  followed  by  a  short  conclusion  in  Section  7.  Subsequent  manuscripts 
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will  describe  the  discretization  and  mapping  of  the  sea  surface  to  the  observer’s  receiver 
optics,  the  geometrical  and  infrared  models  for  the  ship  wake  and  the  performance  of  the 
simulator. 


2  Nomenclature 


Radiation  theory  is  a  mature  field  in  which  has  emerged  a  nomenclature  that  facilitates  a 
precise  description  of  the  physical  quantity  juxtaposed  with  the  interaction  geometry.  In 
Table  1  we  have  described  a  few  of  the  more  commonly  used  radiation  parameters  that  we 
use  in  our  development  of  the  sea  surface  radiance  model.  In  particular  we  will  require  the 
incident  radiation  on  the  sea  surface  from  the  sky,  the  reflected  radiation  and  the  emitted 
radiation  to  model  the  total  radiance  of  the  sea  surface. 


Symbol 

Quantity 

Description 

Units  (SI) 

U 

Radiant  energy 

Energy  traveling  in  the  form  of  EM  waves 

J 

Radiant  flux 

Time  rate  of  flow  of  radiant  energy 

W 

M 

Emittance,  Exitance 

Radiant  flux  emitted  per  unit  area  of  the 
radiating  source 

W  /  nr 

E 

Irradiance 

Total  radiant  flux  incident  on  a  surface  per 
unit  area  of  the  irradiated  surface  per  unit 
of  solid  angle  around  the  irradiative  sur¬ 
face 

W / m 2 

I 

Radiant  Intensity 

Total  radiant  flux  per  unit  of  solid  angle 
around  the  radiating  source 

W/sr 

L 

Radiance,  Stera- 

diancy,  Specific 

intensity,  Brightness 

The  rate  of  energy  radiation  (radiant  flux) 
from  a  source  per  unit  of  solid  angle  and 
per  unit  of  area  of  the  source 

W  /  ( sr.m 2) 

Table  1:  Nomenclature  of  some  radiation  parameters. 


3  Framework  for  the  Sea  Surface  Infrared 
Radiance  Simulator 


Calculating  the  infrared  radiance  received  by  an  observer  from  the  sea  surface  is  a  method¬ 
ical  process  which  in  our  simulator  consists  of  four  modules.  These  describe  the  ambient 
or  background  sea  surface,  the  sea  surface  disturbed  by  the  ship  wake,  the  sky  infrared 
radiance  incident  on  the  sea  surface,  and  the  receiver-source  geometry.  In  Figure  1  we 
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have  illustrated  the  procedure  that  is  implemented  in  our  simulator  to  compute  the  infrared 
radiance  at  each  pixel  of  the  receiver  from  the  appropriate  patch  of  sea  surface.  In  the 
flowchart,  circles  correspond  to  algorithms  while  rectangles  denote  input  or  output  data.  In 
total,  there  are  eight  algorithms  of  which  only  two  will  be  described  in  this  paper. 

The  contribution  from  the  background  sea  surface  is  determined  by  its  roughness  and  tem¬ 
perature  given  a  defined  incident  sky  radiation.  The  roughness  is  largely  determined  by 
the  wind  speed  at  the  sea  surface  and  is  generally  described  as  a  probability  distribution  of 
slopes  of  facets  or  patches  of  the  surface.  The  temperature  will  vary  with  season  and  the  un¬ 
derlying  ocean  currents  that  vary  geographically.  Then  given  the  reflectivity  and  emissivity 
of  the  sea  surface  and  their  inherent  dependencies  on  angle  to  the  observer,  temperature, 
and  roughness  we  can  calculate  the  infrared  radiance  from  the  sea  surface.  These  elements 
are  described  in  detail  in  Sections  4  and  5. 

A  ship  wake  imprints  onto  the  sea  surface  a  distinct  ‘V’  shaped  pattern.  Within  the  ‘V’ 
shape,  the  roughness  of  the  sea  surface  is  sigificantly  modified,  while  outside,  the  back¬ 
ground  sea  surface  persists.  With  distance  downstream  from  the  ship  the  ‘V’  fades  with  the 
background  sea  gradually  infringing  upon  the  wake  until  it  vanishes.  The  wake  consists  of 
the  Kelvin  gravity  waves,  some  white  water  wake  along  the  ship  and  astern,  and  a  trailing 
wake  consisting  of  a  turbulent  bubbly  water  [7].  The  gravity  waves  introduce  a  smooth 
undulation  on  the  sea  surface,  while  the  white  water  is  a  result  of  the  surface  breaking  up 
and  entraining  air.  The  trailing  wake  can  have  roughness  that  is  significantly  different  from 
the  background  sea  and  its  reflectivity  is  thus  different  [8,  9].  Additionally,  if  a  thermal 
stratification  is  disturbed  by  the  ship,  the  trailing  wake  can  be  thermally  different  from  the 
undisturbed  background  sea  surface  [10].  See  Reference  [11]  for  an  introduction  to  various 
elements  of  a  ship  wake  in  the  infrared.  In  our  simulator,  we  will  develop  a  model  for  the 
ship  wake  roughness  and  ship  wake  temperature  which  will  be  published  in  a  subsequent 
paper. 

In  addition  to  the  emitted  radiation  from  the  sea,  infrared  radiation  from  the  sky  which  is 
incident  on  the  sea  surface  is  also  reflected  toward  the  observer.  The  sky,  from  an  observer’s 
point  of  view,  consists  of  the  integrated  effect  of  infrared  emitters  between  the  observer  and 
space.  This  includes  the  celestial  sources  such  as  stars,  the  moon  and  in  particular  the  sun, 
as  well  as  clouds  and  other  emissions  from  atmospheric  molecular  components.  There  are 
several  models  for  skylight  from  clear  skies  [12,  13],  whereas,  cloudy  skies  are  usually 
modeled  with  the  MODTRAN  atmospheric  propagation  model  [14,  15].  In  our  simulator, 
we  expect  to  integrate  existing  models  for  sky  radiation  that  is  incident  on  the  sea  surface. 
In  a  future  paper,  we  wifi  describe  the  choice  of  sky  model  and  its  implementation. 

The  receiver-source  geometry  model  will  accurately  map  a  discrete  facet  of  the  sea  sur¬ 
face  to  a  pixel  in  the  receiving  element.  The  mismatch  between  the  sea  surface  which  is 
an  extended  source  on  a  curved  earth  and  a  rectangular  pixel  at  the  receiver  results  in  a 
non-trivial  discretization  of  the  sea  surface  in  the  numerical  calculations.  These  will  be 
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Figure  2:  The  co-ordinate  system  and  the  geometry  of  a  sea  surface  facet. 

discussed  in  a  subsequent  manuscript.  With  these  four  modules,  we  intend  to  develop  an 
infrared  sea  surface  radiance  simulator  capable  of  modeling  the  infrared  radiance  of  a  ship 
wake. 


4  Sea  Surface  Roughness  and  Reflectivity 


In  modeling  the  sea  surface,  we  consider  it  to  be  tessellated  by  planar  facets.  The  facet 
geometry,  shown  in  Figure  2,  is  the  basic  element  of  the  sea  surface  for  which  we  will 
calculate  the  reflected  and  emitted  radiances.  Planar  facets  will  fully  cover  the  sea  surface 
with  each  facet  approximating  the  orientation  of  the  local  sea  surface  up  to  the  first  order 
derivative  or  slope  in  two  dimensions. 

In  order  to  define  the  sea  surface  model  we  require  a  distribution  of  facet  slopes  that  is 
representative  of  the  real  sea.  There  are  several  experimental  measurements  and  therefore 
experimentally-derived  models  of  the  distribution  of  sea  surface  slopes  under  various  wind- 
driven  sea  states  e.g.  the  Cox  &  Munk  model  [5],  the  Shaw-Churnside  model  [16]  etc.  In 
this  development  we  will  use  the  Cox  &  Munk  (CM)  model,  however,  we  are  able  to 
implement  other  sea  surface  models.  The  CM  model  for  the  probability  of  a  facet  with 
slope  z.x  and  zy  when  the  sea  surface  is  driven  by  a  wind  of  speed  w  m/s  is  given  by 


(1) 
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where 


<52u  =  3.16  x  10~3w  (2) 

<32  =  0.003  +  1.92  x  10'3w 

In  Eq  1,  zx  is  the  slope  of  the  facet  in  upwind  direction  while  zy  is  the  slope  in  crosswind 
direction.  The  mean  value  of  the  slopes  is  zero  while  their  variances  in  the  upwind  and 
crosswind  direcstions  are  (52  and  a2,  respectively. 

To  calculate  the  total  radiance  received  at  the  receiver  from  a  facet,  we  must  account  for  the 
reflected  and  emitted  components.  Below  we  address  the  reflected  part  while  the  emitted 
component  is  described  in  Section  5.  The  reflected  radiance  at  the  receiver  from  a  facet 
depends  on  the  reflectivity  of  the  facet  in  the  direction  of  the  receiver  and  on  the  incident 
radiance  from  the  part  of  the  source  that  can,  by  geometry,  be  reflected  to  the  receiver. 
As  shown  in  Figure  2  the  normal  to  the  facet  and  the  directions  to  the  receiver  and  source 
are  in  standard  spherical  co-ordinates  with  polar  or  zenith  angle  0  and  azimuthal  angle  (|). 
Simple  geometry  demands  that  specular  reflection  off  the  facet  with  normal  (0„.(|)„)  that 
is  received  at  a  receiver  in  the  direction  (0r,  <|>r)  must  originate  from  the  source  at  (0V.  (|)v) 
such  that  the  angle  of  incidence  cos  equals  the  angle  of  reflection  co,-. 

In  our  model,  the  direction  of  the  observer  or  receiver  (0r,  <|)r)  and  the  slopes  of  the  facet 
(zx,zy)  are  known.  The  directions  of  the  facet  normal  (0, ,.())„)  and  the  source  (0y.(|)y)  have 
to  be  expressed  in  terms  of  these  known  parameters.  A  little  algebra  gives 


and 


(3) 


atan 

(  Zy 
Ui, 

zx  <  0  and  zy  <  0 , 

(]),7  —  < 

atan 

(  Zy 
Ui, 

+  2,71 

zx  <  0  and  zy  >  0 , 

(4) 

atan 

(  Zy 
\Zx  j 

+  71 

Zx  >  0. 

In  Figure  3  we  have  plotted  the  orientations  of  the  facet  normal  vector  h  —  (0„.  ())„)  for 
the  range  of  different  slopes.  Equations  3  and  4  express  the  direction  of  the  facet  normal 
(0,7,  (f)/z)  in  terms  of  the  presumed  facet  slopes  z*  and  zy. 
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Figure  3:  The  orientation  of  the  normal  vector  h  =  (0„, <|)„)  for  a  sea  surface  facet  as  a 
function  of  the  slopes  zx  and  zy. 
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The  part  of  the  sky  source  that  is  reflected  toward  the  receiver  is  determined  by  oov  =  co,-  = 
co.  Note  that  since  n.  s  and  r  (the  unit  vectors  for  the  facet  normal,  source  and  receiver)  are 
in  the  same  plane  we  find  that 


co505  =  2  t,cosQn  —  cosQr ,  (5) 

and 

cos§s  —  2C,sinQncos§n  —  sbiQ,  cos§rsinQs ,  (6) 

where 

C,  =  S'//iQ/iCos(j)/l\niQrco\(j)r  +  sinQnsin§nsindrsin§r  +  cosQncosQr  —  n.r  —  cost 0  (7) 


Equations  5-7  express  the  orientation  of  the  vector  to  the  part  of  the  sky  that  is  reflected  by 
a  facet  with  slope  (zx-Zy)  toward  a  receiver  at  (0,-, (]),.).  This  completes  the  basic  geometric 
relations  for  calculating  the  reflected  component. 

We  denote  the  sky  or  source  radiance  incident  on  a  facet  from  the  direction  (0v,c|)s)  to  be 
Ls.  Then  the  reflected  radiance  in  the  direction  of  the  receiver  L,-(0V.(|)V)  is  given  by 

M6rA)  =  P  (0r,4>r)  xL,(0,,4>,)  ,  (8) 

where  the  p(0,-,(|),  )  is  the  reflectivity  (reflectance)  in  the  direction  of  the  receiver.  We 
specify  the  reflectivity  of  the  sea  water  surface  in  Section  6.  Equation  8  is  a  local  relation 
for  the  incident  radiation  reflected  toward  a  receiver  by  a  single  facet.  In  contrast,  a  pixel 
at  the  receiver  will  likely  sample  an  area  of  sea  surface  that  is  composed  of  many  facets. 
Therefore,  the  reflected  radiance  at  the  receiving  pixel  will  be  the  additive  contribution 
of  the  reflection  off  many  facets.  In  fact  we  will  assume  that  the  number  of  facets  is 
sufficiently  numerous  so  that  their  slopes  are  distributed  according  to  the  Cox  &  Munk 
probability  function  of  Eq.  1 .  We  denote  the  horizontal  projection  of  the  area  of  all  facets 
with  slope  (zx ± d Zx-Zy  ± dzy)  by  dA /,  (See  Fig.  2  for  a  schematic  for  a  single  facet).  Adding 
the  area  projections  over  all  the  facets  that  are  visible  to  the  receiver  pixel  we  obtain  the 
total  horizontally  projected  area  A/,.  Then  the  ratio 


dA] i 
Ah 


p(zx,zy)  dzx  dZy, 


(9) 


is  the  fractional  occurence  of  facets  with  slopes  (zx  ±  dzx.  zy  ±  dzy)  that  are  viewed  by  the 
receiver  pixel.  From  this  basic  relation,  it  follows  that 


dA 


n 


Ah 


^p(z„Zy)dz,dzy, 


(10) 
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and 


dAr 

Ah 


COS(£>i- 

COsQn 


P(zX,Zy)  dZxdZy, 


(11) 


which  is  the  same  relationship  but  expressed  in  terms  of  area  projections  on  the  facet  normal 
and  toward  the  receiver.  Denoting  by  Ar  the  total  area  projected  by  the  facets  toward  the 
receiver  we  have 


dAr 

a7 


A/,  cos co, 
Ar  cosQn 


P(zX,Zy )  dzx  dZy. 


(12) 


Equation  8  describes  the  contribution  to  the  reflected  radiance  from  the  sky  to  the  receiver 
from  a  single  facet  with  slope  (zx,zy)  while  Eq.  12  defines  the  area  weight  with  which 
facets  with  slopes  (zx  ±dzX:zy  ±  dzy)  are  projected  toward  the  receiver.  Then  it  follows  that 
the  reflected  radiance  at  the  receiver  pixel  is  given  by 

A.]  cOi^co 

dLr(QrAr)  =  p(0r,<M  x  MOnP.*)  x  — — — ^ p(zx,zy )  dzx  dzy .  (13) 

COS\j pi 


Note  that  when  co,  >  90,  the  facet  is  pointing  away  from  detector.  So,  only  slopes  with 
co,-  <  90  should  be  considered. 


Finaly,  with  po=  j1  p{zx,zy )  dzx  dzy  as  a  normalisation  factor  for  the  probability 


co<! 

distribution  function,  the  total  reflected  radiance  arriving  at  a  single  pixel  at  receiver  is  then 
the  sum  over  all  slopes: 


Lr(Qn$r)  =  — 
P  o 


P  (Or,  tyr)  X  Ls  (0s,  p.y)  X 


A/,  COS  CD, 
A,-  cosQn 


p{zX;Zy)  dzx  dz 


■y- 


(14) 


5  Sea  Surface  Emissivity 


The  derivation  of  the  total  emitted  radiation  from  the  sea  surface  that  is  received  is  similar 
to  that  in  Section  4.  From  a  single  facet,  the  infrared  radiation  emitted  toward  a  receiver  in 
the  direction  (0,-,  <|)r)  is  given  by 


(On  (V) 


£  (On  Pr) 


PjT^) 

n 


(15) 


DRDC  Atlantic  TM  2010-280 


9 


with 


P(T,k) 


2nhc2 

.  c  .  ch  ’ 

X5(e^T  -  l) 


(16) 


where  £  (0,-,  §r)  —  1  —  p  (0r,  §r)  is  the  emissivity  and  P(T.  A.)  is  the  Planck  black  body  func¬ 
tion  at  temperature  T  (c,  h.  and  k  are  physical  constants  for  the  speed  of  light,  Planck’s 
constant  and  Boltzmann  constant). 

A  pixel  at  the  receiver  is  radiated  by  many  facets  with  the  different  slopes.  Once  again 
Eq.  12  quantifies  the  area  weight  projected  toward  the  receiver  of  the  facets  with  slopes 
(zx  ±  dz.x-Zy  ±  dzy)  .  It  follows  that  the  emitted  radiation  received  from  these  facets  is 


dLe{Qr,$r)  =  [1  —  p (0r,4>r)]  x  x  ^lcos(^'  p(zx,Zy)  dzx  dzy .  (17) 

7C  r\.f  COSyjfi 

The  total  emitted  radiance  at  a  receiver  pixel  is  given  by  integrating  over  all  slopes: 


,n  ,  .  1  f  f  ..  P(T,X)  Ah  cosar  .  . 

Le(Qr,$r)  =  —  /  /  [1  —  p  (0r,  (>r)]  X - -  x  - —p{zx,Zy )  dzx  dzy . 

PO  v  ^  ^  ^ r  COS\Jn 

G><! 


(18) 


6  Results  and  Discussion 


Our  principal  result  is  the  total  radiance  at  a  receiver  pixel  from  a  patch  of  sea  surface  that 
is  comprised  of  many  facets.  The  total  radiance  is  the  sum  of  the  reflected  (Eq.  14)  and 
emitted  (Eq.  18)  contributions: 


L(0,.,4v)  =  Le  (0,.,  §r)  +  L,  (0,.,  (]),.) 


1 

Po 


[1  —  p  (0r,<>r)]  x 


P(T,  X) 


n 


co<5 


T  P  (0rj  tyr)  x  Ls  (0J,  (fh)}  x 


A/t  cos(Or 
Ar  cosQn 


P{Zx,Zy ) 


dzx  dZy 


(19) 
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This  radiance  model  for  the  Cox  &  Munk  sea  surface  is  comparable  to  several  others  that 
are  currently  used  in  the  scientific  community,  see  Fig.  4  for  a  short  list. 

In  calculating  Eq.  19,  we  map  a  patch  of  sea  surface  to  a  receiver  pixel  given  the  instan¬ 
taneous  field  of  view  of  the  pixel  and  the  viewing  position  and  orientation.  The  patch  of 
sea  surface  is  covered  by  a  collection  of  facets  whose  slopes  are  sampled  from  the  Cox 
&  Munk  probability  distribution.  For  each  patch  we  determine  the  orientation  toward  the 
receiver  and  then  compute  the  emitted  and  reflected  radiation.  For  these  we  need  to  specify 
the  emissivity  or  reflectivity  p(0, ,  <|)r)  for  the  sea  water  facet  and  the  incident  sky  radiation 
LV(0V,  (|)v).  The  reflectance  can  be  determined  from  the  index  of  refraction  of  water  by  the 
Fresnel  equations.  In  Table  2  we  have  listed  the  index  of  refraction  of  seawater  for  infrared 
radiation  of  wavelength  3.4  —  5.2/rm  at  25°C  from  Hale  &  Querry  [17].  The  resulting  re¬ 
flectivity  and  emissivity  are  independent  of  <[>,-  and  vary  with  the  zenith  angle  0;.  as  shown 
in  Fig.  5.  When  viewed  from  directly  overhead  (0,  =  0),  the  sea  surface  is  essentially  emis¬ 
sive  with  near- vanishing  reflectance,  while  at  grazing  angles  (0r  ~  90°)  the  sea  surface  is 
primarily  reflective  with  near-vanishing  emissivity.  The  strong  variation  of  the  reflectance 
with  zenith  angle  is  expected  to  significantly  impact  the  total  radiance  at  the  receiver. 


X  [m] 

nfA,) 

X  [m] 

nfA.) 

X 

[m] 

nfA.) 

X 

[m] 

n(A) 

3.4 

1.42 

3.9 

1.357 

4.4 

1.334 

4.9 

1.328 

3.5 

1.4 

4 

1.351 

4.5 

1.332 

5 

1.325 

3.6 

1.385 

4.1 

1.346 

4.6 

1.330 

5.1 

1.322 

3.7 

1.374 

4.2 

1.342 

4.7 

1.330 

5.2 

1.317 

3.8 

1.364 

4.3 

1.338 

4.8 

1.330 

5.3 

1.312 

Table  2:  The  index  of  refraction  of  water  for  infrared  light  at  25°  C. 


The  sky  radiance  that  is  incident  on  a  facet  of  the  sea  surface  originates  primarily  from 
scattered  sunlight  and  thermal  emission  from  clouds  and  the  atmosphere.  Other  stellar 
sources  such  as  moon  light  are  important  on  clear  nights.  Whereas  the  sky  radiance  is 
heterogeneous,  for  the  purposes  of  demonstrating  our  simulator,  we  have  chosen  a  sky 
radiance  that  is  independent  of  the  azimuthal  and  zenith  angles.  We  illustrate  two  cases 
where  the  constant  sky  radiance  is  greater  or  less  than  the  thermal  emission  from  the  sea 
surface.  Given  the  sky  radiance  and  the  reflectance  function,  we  can  calculate  the  total 
radiance  at  the  receiver  for  a  given  sea  temperature  and  sea  surface  roughness. 

In  Fig.  6  we  have  plotted  the  emitted,  reflected  and  total  sea  surface  radiance  at  the  reciver 
for  a  sky  radiance  that  is  greater  and  lesser  than  the  emitted  sea  radiance.  We  find  that  the 
total  radiance  increases  with  the  zenith  angle  generally  following  the  shape  of  the  reflected 
sea  surface  radiance  which  is  increasingly  important  at  grazing  angles.  When  the  sky 
radiance  is  lower  than  the  intrinsic  sea  surface  radiance,  the  total  radiance  decreases  with 
the  zenith  angle  generally  following  the  shape  of  the  intrinsic  sea  surface  radiance.  As 
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noted  earlier,  these  behaviours  stem  from  the  dependence  of  the  emissivity  and  reflectance 
on  the  zenith  angle  as  determined  by  the  Fresnel  equations. 

At  constant  temperature,  the  variation  of  the  reflected  radiance  with  sea  roughness  deter¬ 
mines  the  variation  of  the  total  sea  radiance.  In  Fig.  7  we  have  plotted  the  total  radiance 
received  as  a  function  of  wind  speed,  effectively  sea  surface  roughness,  for  a  zenith  an¬ 
gle  of  75°  for  an  incident  sky  radiance  that  exceeds  the  emitted  sea  surface  radiance.  At 
grazing  angles,  the  angle  between  the  receiver  and  the  normal  to  the  facets  decreases  with 
increasing  roughness.  This  is  equivalent  to,  on  average,  a  lower  reflectvity  and  a  higher 
emissivity.  Consequently,  when  the  sky  radiance  is  greater  (lesser)  than  the  emitted  sea 
surface  radiance,  the  total  sea  surface  radiance  increases  (decreases)  with  roughness.  At 
non-grazing  angles,  the  angle  between  the  receiver  and  the  normal  to  the  facets  increase 
with  roughness.  On  average,  this  amounts  to  an  increase  in  reflectivity  and  a  decrease  in 
emissivity.  Consequently,  when  the  sky  radiance  is  greater  (lesser)  than  the  emitted  sea  sur¬ 
face  radiance,  the  total  sea  surface  radiance  increases  (decreases)  with  roughness  as  shown 
in  Fig.  8. 

In  contrast  to  the  variation  of  sea  surface  radiance  with  roughness,  the  dependence  of  the 
emitted  radiation  on  sea  surface  temperature  is  monotonic:  in  Fig.  9,  we  have  illustrated 
this  behaviour  for  grazing  and  non-grazing  angles  and  for  sky  radiances  greater  than  and 
lesser  than  the  sea  surface  radiance. 

7  Conclusion  and  Future  Work 


In  this  manuscript  we  have  described  the  framework  and  development  of  a  sea  surface  radi¬ 
ance  simulator  that  is  envisioned  to  account  for  the  effects  of  optically  coherent  structures, 
particularly  ship  wakes  on  the  sea  surface,  as  well  as  the  pseudo-random  wind  driven  sea 
surface.  Our  approach  is  designed  to  utilize  existing  third  party  modules  for  modeling 
the  radiance  from  the  sky,  atmosphere  and  propagation  of  infrared  radiation.  At  this  time 
we  have  described  the  model  for  the  background  sea  surface.  It  consists  of  a  wind  driven 
surface  that  can  be  locally  approximated  by  planar  facets  defining  spatial  gradients  in  the 
upwind  and  cross  wind  directions.  The  orientation  of  the  facets  or  more  directly  the  slopes 
that  determine  the  orientation  are  sampled  from  an  empirical  probability  distribution.  With 
well  defined  reflectance  and  emissivity  for  sea  water,  we  modeled  the  reflected  sky  radi¬ 
ance  and  the  emitted  black  body  radiance  by  the  sea  surface.  With  our  simulation  we  have 
illustrated  the  dependence  of  the  sea  surface  radiance  on  roughness  and  temperature  which 
serve  to  verify  the  functioning  of  the  background  sea  surface  simulator. 

Subsequent  manuscripts  will  describe  three  other  modules  that  interface  with  the  back¬ 
ground  module  described  here.  A  sky  radiance  model  based  on  third  party  software  ShipIR 
and  atmospheric  propagation  software  MODTRAN  will  be  incorporated  to  take  into  ac¬ 
count  realistic  distributions  of  sky  radiance  from  solar  scattering,  non-uniform  cloud  cover, 
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aerosols  and  stellar  radiation.  A  home-grown  model  mapping  the  extended  sea  surface  to 
a  receiver  will  be  published.  The  model  describes  the  receiver  optics  relating  the  signal 
received  at  a  single  pixel  to  the  patch  of  sea  surface  that  is  imaged,  and  the  subsequent 
discretization  of  the  patch  into  a  collection  of  planar  facets.  And  finally,  a  module  for  the 
ship  wake  is  under  development.  The  model  will  account  for  the  various  features  of  a  ship 
wake  taking  into  account  the  undulation  of  the  sea  surface  due  to  the  Kelvin  gravity  waves, 
the  damped-capillary  waves  in  the  trailing  wake,  the  thermal  wake  in  stratified  waters  and 
the  white  water  wake.  The  capability  to  accurately  model,  in  the  optical  and  infrared  bands, 
a  ship  wake  in  the  background  sea  is  at  present  a  significant  deficiency  in  the  infrared  sea 
and  target  modeling  community. 
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Figure  4:  Models  for  the  sea  surface  roughness  and  radiance. 
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Reflactance  and  Emissivity 


Figure  5:  The  dependence  of  emissivity  and  reflectivity  of  the  sea  surface  on  the  angle  of 
incidence. 
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Sea-surface  radiance  depending  on  the  Zenith  angle  for  Tsea=6.5C 


Figure  6:  The  emitted,  reflected  and  total  sea  surface  radiance  as  a  function  of  the  zenith 
angle  for  a  uniform  sky  with  sky  radiance  greater  than  the  emitted  sea  surface  radiance 
(top)  and  sky  radiance  less  than  the  emitted  sea  surface  radiance  (bottom). 


16 


DRDC  Atlantic  TM  2010-280 


a) 


Roughness  (wind  speed)  [m/s] 


Figure  7:  The  dependence  of  the  sea  surface  radiance  on  sea  surface  roughness  for  a  zenith 
angle  of  75°  with  (a)  constant  sky  radiance  greater  than  the  emitted  sea  radiance  and  (b) 
constant  sky  radiance  lesser  than  the  emitted  sea  radiance. 


DRDC  Atlantic  TM  2010-280 


17 


a) 


Figure  8:  The  dependence  of  the  sea  surface  radiance  on  sea  surface  roughness  for  a  zenith 
angle  of  25°  with  (a)  constant  sky  radiance  greater  than  the  emitted  sea  radiance  and  (b) 
constant  sky  radiance  lesser  than  the  emitted  sea  radiance. 
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al) 


Figure  9:  The  dependence  of  the  sea  surface  radiance  on  the  water  temperature  for  a  uni¬ 
form  sky  with:  al)  sky  radiance  greater  than  the  emitted  sea  surface  radiance  for  a  high 
zenith  angle;  bl)  for  a  low  zenith  angle;  a2)  sky  radiance  lesser  than  the  intrinsic  sea  sur¬ 
face  radiance  for  a  high  zenith  angle;  b2)  for  a  low  zenith  angle. 
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